Quark models of dibaryon resonances in nucleon-nucleon scattering 



J.L. Ping 1 , H.X. Huang 1 , H.R. Pang 2 , Fan Wang 3 and C.W. Wong 4 
1 Department of Physics, Nanjing Normal University, Nanjing, 210097, P.R. China 
2 Department of Physics, Southeast University, Nanjing, 210094, P.R. China 
3 Department of Physics, Nanjing University, Nanjing, 210093, P.R. China 
4 Department of Physics and Astronomy, University of California, Los Angeles, CA 90095-1547, USA 

We look for AA and NA resonances by calculating AW scattering phase shifts of two interacting 
baryon clusters of quarks with explicit coupling to these dibaryon channels. Two phenomenological 
nonrelativistic chiral quark models giving similar low-energy NN properties are found to give sig- 
nificantly different dibaryon resonance structures. In the chiral quark model (ChQM), the dibaryon 
system does not resonate in the NN S-waves, in agreement with the experimental SP07 NN partial- 
wave scattering amplitudes. In the quark delocalization and color screening model (QDCSM), the 
5*-wave NN resonances disappear when the nucleon size b falls below 0.53 fm. Both quark models 
give an IJ P = 03 + AA resonance. At b — 0.52 fm, the value favored by baryon spectrum, the 
resonance mass is 2390 (2420) MeV for the ChQM with quadratic (linear) confinement, and 2360 
MeV for the QDCSM. Accessible from the 3 Df N channel, this resonance is a promising candidate 
for the known isoscalar ABC structure seen more clearly in the pn^d-Kix production cross section 
at 2410 MeV in the recent preliminary data reported by the CELSIUS- WASA Collaboration. In 
the isovector dibaryon sector, our quark models give a bound or almost bound 5 S^ A state that can 
give rise to a 1 D2 N resonance. None of the quark models used has bound NA P-states that might 
generate odd-parity resonances. 



PACS numbers: 12.39.Jh, 14.20.Pt, 13.75.Cs 

I. INTRODUCTION 

Quantum chromodynamics (QCD) is widely accepted 
as the fundamental theory of the strong interaction. Lat- 
tice QCD methods have recently been used to study 
low-energy hadronic interactions, including the nucleon- 
nucleon (NN) interaction [l[. However, QCD-inspired 
quark models are still the main tool for detailed studies 
of the baryon-baryon interaction. 

The phenomenological quark model most commonly 
used in the study of NN interaction is the nonrelativis- 
tic chiral quark model (ChQM) @, 0, H- Nonrelativis- 
tic kinematics makes the many-body treatment of the 
multiquark system manageable, with the very convenient 
choice that the light quark mass m q is just one third of 
the nucleon mass. With quarks, one needs a confinement 
potential to reproduce a distinctive QCD property, and 
a one gluon exchange (OGE) to account for A—N mass 
difference and other details of baryon excitations. The 
inclusion of pion exchange will take care of long-range 
baryon-baryon interactions as well as some features in 
baryon structure, and is the consequence of the relative 
weakness of chiral symmetry breaking. Finally, scalar 
exchange is used to describe an extra intermediate-range 
attraction needed in nuclear forces. No other meson ex- 
changes are included. ChQMs with a few chosen and a 
few adjusted parameters are able to give a surprisingly 
simple if only semi-quantitative picture of both baryon 
spectra and baryon-baryon interactions at relatively low 
energies. It is therefore of some interest to understand 
some of the limitations of these simple quarks models. 

The most problematic term in the ChQM is the scalar 
exchange term. It takes into account neglected chan- 



nels containing As and pions, and should therefore vary 
when more of these channels are explicitly included in 
the coupled-channel calculation. Its effects always in- 
clude the exchange of two pions, and are called correlated 
if the two pions also interact with each other. Modern 
treatments of correlated two-pion exchange show that in 
addition to a long-range scalar-isoscalar attraction tradi- 
tionally associated with scalar exchange, there is also a 
strong scalar-isoscalar repulsive core |5|, a complication 
that has not yet been included in the ChQM. 

In the quark model, the forces between baryonic clus- 
ters of quarks and antiquarks are like molecular forces 
between molecules of atoms built up from the forces be- 
tween their constituents. This molecular model of nuclear 
forces has been extensively developed by using the quark 
delocalization and color screening model (QDCSM) [f|. 
In this model, quarks confined in one baryon are allowed 
to delocalize to a nearby baryon and to change the dy- 
namics of the baryon-baryon interaction through a reduc- 
tion of the confinement potential called color screening. 
The delocalization parameter that appears is determined 
by the dynamics of the interacting quark system, thus 
making it possible for the quark system to reach a more 
favorable configuration through its own dynamics in a 
larger Hilbert space. The model has been successfully 
applied to NN and hyperon-nucleon scatterings. The 
important intermediate-range NN attraction is achieved 
in this model by the mutual distortion of the interact- 
ing nucleons, in a way that is very similar to the mutual 
distortion of interacting molecules. 

The main difference between the ChQM and the QD- 
CSM is the mechanism for intermediate-range attraction. 
Recently, we showed that both the ChQM containing the 
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cr-meson and the QDCSM without it gave a good descrip- 
tion of the low-energy A A S- and D-wave scattering 
phase shifts and the properties of deuteron with almost 
the same quark- model parameters Q ■ Thus the a- meson 
exchange can effectively be replaced by the quark der- 
ealization and color screening mechanism. It is not clear 
however if their equivalence persists to higher energies 
where nucleons overlap more strongly and baryon exci- 
tation and multiquark effects become more important. 

Interest in multiquark system has persisted since R. 
Jaffe predicted the H-particle in 1977 .Si]. All quark mod- 
els, including those using lattice QCD techniques, predict 
that in addition to the qq mesons and q 3 baryons, there 
should be multiquark systems (qq)\ q 4 q, q 6 , quark gluon 
hybrids qqg, q 3 g, and glueballs Up to now there 
has been no well established experimental candidate of 
these multiquark states [ltj. Recently, the CELSIUS- 
WASA Collaboration has reported preliminary results 
on the ABC anomaly in the production cross section of 
the pn—>dir°7r° reaction that suggests the presence of an 
isoscalar J p = 1 + or 3 + subthreshold AA resonance, 
with resonance mass estimated at ~ 2410 MeV and a 



width of < 100 MeV 11]. The relatively large binding en- 
ergy involved gives an object that is much closer to these 
interesting multiquark states than a loosely bound sys- 
tem like the deuteron. Nonrelativistic quark models such 
as ChQMs and QDCSMs fitting both A, A masses and 
low-energy A A scattering properties can be expected to 
give particularly interesting and parameter-free predic- 
tions for such dibaryon multiquark states. 

It thus appears worthwhile to extend our past calcu- 
lation of NN phase shifts to the resonance region near 
the AA and A A thresholds by including these excited 
dibaryon channels in coupled-channel calculations. The 
A resonance is by far the most important low-energy 
baryon resonance. It dominates even the ir~p cross sec- 
tions where its production is hindered relative to the pro- 
duction of isospin 1/2 N* resonances by a factor of 2 
from isospin coupling. In dibaryon channels, the AA 
threshold at 2460 MeV is clearly separated from the 
AA*(1440) threshold at 2380 MeV, the second most in- 
teresting dibaryon threshold in this energy region. A A 
bound states are of immediate interest in understanding 
resonance phenomena in this near subthreshold energy 
region. The inclusion of V*(1440) would be of consider- 
able interest in a broader study of dibaryon resonances, 
but its inclusion is technically difficult for us because 
TV* (1440) is commonly understood to be a monopole ex- 
citation of the nucleon that has a much more compli- 
cated quark wave function. In contrast, our approximate 
description of the N — > A excitation as a simple spin- 
isospin excitation without any change in the radial wave 
function probably captures the essence of the physics in- 
volved. A excitations are thus within easy reach of the 
technology used in our previous coupled-channel calcula- 



tions. Our first study in the resonance region will include 
only A excitations. This limitation of excited baryon de- 
grees of freedom to As only has often been made in past 
studies of nuclear forces. [12| 

We shall use the same Salamanca ChQM fl3j and QD- 
CSM used previously for NN channels only Q with addi- 
tional sets of potential parameters to find out if their sim- 
ilarity persists into the resonance region. We are inter- 
ested in particular in discovering how similar these sim- 
ple quark models are in describing theoretical dibaryon 
resonances originating from AA or AA bound states 
when the As are treated as stable particles. In other 
words, these dibaryon resonances are theoretical com- 
pound dibaryon states that are allowed to decay only via 
the NN channels. A brief description of these two quark 
models of the baryon-baryon interaction is given in Sec- 
tion II. 

The NN phase shifts for these quark-quark potentials 
are calculated using a coupled-channel resonating group 
formalism that includes A A, A A and sometimes 
hidden-color channels as well. No explicit pionic channels 
are included as the As are treated as stable particles. In 
this context, a promising dibaryon resonance is taken to 
be one arising from a bound state below the AA or AA 
threshold. The calculated results, including resonance 
masses and widths (FWHMs), are given in section III in 
those partial waves (PWs) where a theoretical dibaryon 
resonance appears in at least one quark model. No AA 
bound state is found in any isovector odd-parity state in 
all our quark models. 

In Sect. IV, these results are compared to partial-wave 
analyses of A A scattering amplitudes [HI, [l6| , where the 
presence of a dibaryon resonance causes a rapid counter- 
clockwise motion in the Argand diagram. The possibil- 
ity that the ABC effect in the pn^xl-KTr reaction is an 
isoscalar A A resonance is also discussed. 

Section V contains brief concluding remarks on what 
we have learned about quark dynamics in the AA reso- 
nance region. 



II. TWO QUARK MODELS OF 
BARYON-BARYON INTERACTIONS 

A. Chiral quark model 

The Salamanca ChQM is representative of chiral quark 
models. It has also been used to describe both hadron 
spectroscopy and nucleon-nucleon interactions. The 
model details can be found in [3]. Only the Hamilto- 
nian and parameters are given here. 

The ChQM Hamiltonian in the baryon-baryon sector 
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Here S^ is quark tensor operator, Y(x), H(x) and G(x) 
are standard Yukawa functions [1 71 ] , T c is the kinetic 
energy of the center of mass, a c h is the chiral coupling 
constant, determined as usual from the 7r-nucleon cou- 
pling constant. An additional p meson exchange poten- 
tial V between quarks has been added to give an im- 
proved treatment of baryon-baryon P-states. Its param- 
eters will be specified in Sect. IIIIl All other symbols have 
their usual meanings. 

Table Q] gives the model parameters used. For each 
set of parameters, the nucleon size b that appears in 
Eq.([5]) is given a pre-determined value. Two of the 
parameters (a c ,a s ) are fitted to A— N mass difference 
(1232 — 939 MeV) and the equilibrium condition for the 
nucleon mass at the chosen b. The absolute nucleon mass 
is controlled by a constant term Vo in the confinement 
potential that does not affect the baryon-baryon interac- 
tion. In ChQM2 [Hj], the deuteron binding energy (2.22 
MeV) is fitted by varying the combination m a , A calcu- 
lated with the standard two NN coupled channels 3 Si 
and 3 Di (called 2NNcc in the following). The remaining 
parameters m q , m^ and a c u are fixed at chosen values. 
The numeral 2 in the name ChQM2 refers to its quadratic 



confinement potential. The model ChQMl Q uses a lin- 
ear confinement potential instead. The model ChQM2a 
differs from ChQM2 in fitting a different equilibrium nu- 
cleon size b but for simplicity m a , A are allowed to re- 
main at the ChQM2 values. Its deuteron binding energy 
is reduced, but we do not consider the difference to be 
important in our study of dibaryon resonance properties. 



Finally, Table U also gives the effective-range (ER) pa- 
rameters of a 5-parameter ER formula in the NN 3 Si 
( 1 5o) state. The calculation uses with 5 (4) color-singlet 
channels denoted 5cc (4cc) and defined in the following 
section. The channels used include 2 NN and 3 AA (an 
NN, 2 AA, and an NA) channels. The deuteron bind- 
ing energy Ed calculated from the two ER parameters of 
the 5cc calculation is also given in the table. Accord- 
ing to , this approximation overestimates the binding 
energy, but by only 0.015 MeV. So the tabulated bind- 
ing energies are sufficiently accurate for this qualitative 
study. The same ER aproximation for a 2NNcc calcula- 
tion gives Ed ~ 1.86 MeV for ChQM2, thus showing that 
the 3 A A channels increase Ed by about 1.5 MeV. 
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TABLE I: Parameters that differ in different models are given 
in this table. The dimension of each dimensional parame- 
ter is given within parentheses following each symbol: b (fm) , 
a c (MeVfm -2 if quadratic, but MeVfm" 1 if linear), V (fm 2 ) 
and /i(fm~ 2 ). Parameters having the same value for all the 
quark models discussed in this paper are m q = 313 MeV, 
m n = 138MeV, A/hc = 4.2fm~\ and m CT = 675 (MeV) 
for ChMQs. The scattering length and effective range calcu- 
lated for each potential are also given: at,rt for the triplet 
state 3 Si and a s ,r s for the singlet state 1 So, all in fm. 
The deuteron binding energy Sd (MeV) is calculated from the 
triplet effective-range parameters. 



ChQM2 ChQM2a QDCSMO QDCSM1 QDCSM3 
(ChQMl) 
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0.518 


0.60 


0.48 


0.518 


0.60 


a c 


46.938 


12.39 
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18.55 




(67.0) 
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-1.297 


0.255 


-1.299 


-1.3598 


-0.5279 


M 






0.30 


0.45 


1.00 


a s 


0.485 


0.9955 


0.3016 


0.485 


0.9955 


a c h 


0.027 


0.027 


0.027 


0.027 


0.027 




(0.0269) 










a t 


4.52 


20.8 


34.9 


5.94 


6.03 


n 


1.56 


2.24 


2.27 


1.75 


1.67 


Ed 


3.35 


0.11 


0.04 


1.75 


1.64 


a s 


-170 


-2.48 


-2.32 


-6.90 


-5.41 


r s 


2.17 


5.42 


4.48 


2.63 


3.56 



B. Quark derealization, color screening model 

The model and its extension were discussed in detail in 
0, H3] . Its Hamiltonian has the same form as Eq. (JTJ) , but 
used with V = and a different confinement potential 

V° ON (rij ) = -a c Xi ■ Xj [/y (r l3 )+V ], (3) 

where /y(ry) = rjj if quarks each in the Gaussian 
single-quark wave function 4> a of Eq. © , are on the same 
side of the dibaryon, i.e., both centered at S/2 or at 
— S/2. If quarks i, j are on opposite sides of the dibaryon, 

/«(»•«) = i(l-e-"^). (4) 

Quark delocalization in QDCSM is realized by assum- 
ing that the single particle orbital wave function of QD- 
CSM as a linear combination of left and right Gaussians, 
the single particle orbital wave functions of the ordinary 
quark cluster model, 

lMS,e) = (&»(S)+e0a(-S))/tf(e), 
M-S,e) = (M-S)+eMS))/iV(e), 
N(e) = Vl + e 2 + 2ee- s2 / 4b2 ; 

«S, - (_L) S/4 e-*<-^ 



Quak delocalization with color screening is an approxi- 
mate way of including hidden-color (h.c.) effects. 

The color screening constant \x in Eq.Q is determined 
by fitting deuteron properties. The parameters of the 
QDCSMi {%= 1,3) used here are those of Set i of [7|, and 
are given again in Table HI For QDCSMO, \i is an esti- 
mated value that has not been fine tuned to the deuteron 
binding energy. These models differ in the equilibrium 
nucleon size b. 



III. RESULTS 

NN scattering phase shifts are calculated for the quark 
models of Table fl] to energies beyond the AA or iVA 
threshold for different choice of coupled channels. We in- 
clude channels containing one or more As treated as sta- 
ble particles, and channels containing hidden-color (h.c.) 
states. The resonating-group method (RGM), described 
in more details in is used. 

Past experience has suggested that reliable estimates 
of resonance masses can be made using non-decaying As 
0. For the theoretical NA resonance d! (IJ P = 20") 
at the theoretical mass 2065 MeV (and an NA binding 
energy of 106 MeV), an increase in the imaginary part 
of the A resonance energy by 10 — 15 MeV is known 
to increase the d! mass by only a few tenth of an MeV 
pll | . This result suggests that the complete neglect of the 
imaginary part T/2 « 60 MeV of the A resonance energy 
will underestimate the d! resonance mass by perhaps 2 
MeV. If this result holds generally for other resonances 
at other binding energies, our estimates of the resonance 
masses can be expected to be good to a few MeV. 

The use of coupled channels containing stable As does 
mean that the calculated NN phase shifts do not de- 
scribe inelasticities correctly. Thus they cannot be com- 
pared quantitatively to experimental phase parameters 
above the pion-production threshold in NN channels 
with strong inelasticities. For this reason, the primary 
emphasis of this paper is the extraction of resonance en- 
ergies from phase shifts in the resonance region. 

Our theoretical dibaryons are made up of two stable 
constituents below their breakup threshold and are there- 
fore real resonances in the model. They have finite widths 
that come from the coupling to open NN channels. 

A. 1=0 states 

Calculational details and results for the IJ P = 03 + 
states are given in Table |ITJ The number of chan- 
nels used in the theory is given by N c h- The theo- 
retical pure 7 binding energy is next estimated by 
diagonalizing the Hamiltonian matrix for this state in 
a generator-coordinate representation where the average 
baryon-baryon separation is taken to be less than 6 fm (in 
order to keep the matrix dimension manageably small). 
In this way, the pure 7 S^ A is found to be bound by 
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TABLE II: The AA or resonance mass M and decay width Y, 
in MeV, in five quark models for the IJ P = 03 + states. The 
channels included are one channel or lc ( 7 S AA only), two 
coupled channels or 2cc (lc + 3 D$ N ), 4cc (2cc + 7 ' 3 D AA ), 
and 10 coupled channels as described in the text. The pure 
7 S AA bound state mass for ChQMl is 2456 MeV 0. 

~N7 h ChQM2 ChQM2a QDCSMO QDCSM1 QDCSM3 
M Y M Y M Y M Y M Y 

lc 2425 - 2430 - 2443 - 2365 - 2276 

2cc 2428 17 2433 10 2416 20 2368 20 2278 19 

4cc 2413 14 2424 9 2400 14 2357 14 2273 17 

lOcc 2393 14 ______ 

lOcc' 2353 17 ______ 

lOcc" 2351 21 ______ 



35-190 MeV, with the ChQM2 mass 60 MeV lower than 
the ChQMl mass obtained by @. Coupling to the 3 Df N 
channel causes this bound state to change into an elastic 
resonance where the phase shift, shown in Fig. [IJ rises 
through 7r/2 at a resonance mass that has been shifted 
up by 3 MeV. The same small mass increase caused by 
the coupling to the N N continuum is seen in all J p states 
studied here. The result shows that the mass shift is al- 
ways dominated by the NN scattering states below the 
pure bound state mass rather than those above it. 

The table next shows that on coupling to the two 
7 ' 3 D AA channels above the pure 7 S AA bound state, the 
resonance is pushed down in mass, as expected. The ef- 
fect is not large, however, being 15 MeV in ChQM2 and 
11 MeV in QDCSM1, which has the same a s . 

Both the pure 7 S AA bound state and the associated 
3 D^ N resonance are lower in mass in QDCSM1, by about 
60 MeV, than in ChQM2. Since the QDCSMs contain 
h.c. effects, we must first determine how much h.c chan- 
nels contribute to this mass difference. This study is 
done with ten channels (or case lOcc). Besides the four 
baryon channels shown in Table [Til they include the fol- 
lowing six h.c. channels: four 3 __>3 channels of 2 A8 2 As, 
4 N 8 4 N 8 , 4 N 8 2 N 8 , 2 N 8 2 N 8 , the 7 S 3 { 4 N 8 4 N 8 ) and the 
7 D 3 ( 4 N 8 4 N 8 ). Here the baryon symbol is used only to 
denote the isospin so that 2 Ag means the T,S = 3/2, 1/2 
color-octet state. The table shows that these six h.c. 
channels lower the ChQM2 resonance mass by 20 MeV. 
Assuming that the QDCSMs account adequately for h.c. 
contributions, we see that the QDCSM1 dibaryon reso- 
nance mass is now lower than that in ChQM2 by only 36 
MeV. 

It is interesting that the pure 7 S AA bound state masses 
in the ChQM2 and ChQM2a of different nucleon size b 
differ by only 5 MeV. In contrast, both bound state and 
resonance masses change by 90 MeV in the two QDCSMs 
with the same difference in the nucleon size b. This shows 
that the QDCS mechanism depends sensitively on the 
nucleon size. The possibility of a 7 S AA bound state giv- 
ing rise to a subthreshold resonance has been suggested 
previously 22], but the resonance masses for the present 
quark models are higher. 



TABLE III: The AA or resonance mass M and decay width Y, 
in MeV, in five quark models for the IJ P = 01 + states. The 
channels included are lc ( 3 S A A only), 2cc (lc + 3 Sf JV ), and 
5cc (2cc + 3 D? N + 3 - 7 D AA ). The pure 3 Sf A bound-state 
mass for ChQMl is 2274 MeV 0|. 



N ch ChQM2 ChQM2a QDCSMO QDCSM1 QDCSM3 





M 


M 


M 


M 


M 


r 


lc 


2366 


2344 


2317 


2206 


2115 




2cc 


nr" 


nr 


nr 


nr 


2408 


74 


5cc 


nr 


nr 


nr 


nr 


2393 


70 



"No resonance in these coupled channels. 

The calculated 3 D 3 N phase shifts are shown in Fig. Q] 
They differ noticeably between two quark models of very 
similar parameters that differ only in the replacement of 
the scalar exchange in ChQM2 by the QDCS mechanism 
in QDCSM1. The ChQM2 phase shifts agree better with 
experiment than the QDCSM1 phase shifts for 100 < 
E cm < 400 MeV. The situation is different in the 3 S"f N 
state shown in Fig. [5] where the phase shifts from these 
two quark models agree, as already pointed out by 0]. 

Dibaryon resonance masses can be quite sensitive to 
short-range dynamics. For example, the ChQM2 reso- 
nance mass can be forced down to the lower QDCSM1 
value by artificially increasing the attractive confinement 
interaction involving h.c. configurations that appear only 
for overlapping baryons. The change can be made in 
two different ways: Either (1) increase the color confine- 
ment interaction strength among h.c. channels, and also 
that between color-singlet channels and h.c. channels 
(a model denoted lOcc') by an overall factor of 1.3; or 
(2) increase the color confinement interaction between 
color-singlet channels and h.c. channels only (model 
lOcc") by an overall factor 1.55. Fig. \V[b) shows that 
the resulting 3 D 3 N phase shifts change noticeably only 
for E cm > 350 MeV. In other words, only these "high- 
energy" NN phase shifts are sensitive to changes in 
short-range dynamics when shielded by the strong cen- 
trifugal barrier in the NN D-waves. Hence these changes 
in the confinement interactions, artificial as they are, are 
not excluded by the experimental phase shifts. 

The resonance widths given in Table ITT1 are FWHMs. 
They are quite small, and agree with one another. 

Results for the IJ P = 01 + state are shown in Ta- 
ble HITJ and Fig. H The theoretical pure 3 S"f A state is 
bound by 100 - 350 MeV, around twice the 7 S 3 AA bind- 
ing energy. The coupling to the 3 S^ N channel has an 
unexpectedly large effect, pushing up the lowest of these 
bound 3 S AA masses, by 293 MeV in QDCSM3, so that it 
becomes a resonance at 2408 MeV. This very large mass 
shift is caused by the presence of a lower-mass state, 
the deuteron, in the admixed 3 S^ N channel. Admix- 
ing three additional channels with no lower bound states 
pushes the resonance mass down a little to give for the 
5cc treatment the mass shift 

AM ee M R (5cc) - M(lc) = 278 MeV. (6) 
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FIG. 1: (a) 3 D$ N phase shifts calculated for a single channel 
(sc) and two coupled channels (2cc) ( 3 D^ N + 7 S AA ) in two 
different quark models (ChQM2 and QDCSMI ) as functions 
of the cm. kinetic energy E cul = W — 2Mn, where W is 
the total cm. energy and Mn is the nucleon mass. The 
experimental phase shifts of the partial-wave solution SP07 
are also shown, (b) NN 3 D3 phase shifts calculated in 
the ChQM2 using different numbers of coupled channels, as 
described in more detail in the text. 



The pure 3 S AA bound state mass appears 100 MeV or 
more higher in the other four quark models. The addi- 
tional large mass increase caused by the coupling to the 
3gNN channgi without or with the additional channels 
then pushes the state above the AA threshold. Then no 
resonance appears. 

The pure 3 S AA bound state mass in ChQMl is higher 
by 160 MeV than the QDCSM3 mass. So the 3 Sf N 
resonance also does not appear in ChQMl. We shall 
find that QDCSM3 model has an unusually rich dibaryon 
spectrum arising from an unusually strong attraction in 
AA channels. 

It is clear from Table Hill that the strong AA attraction 
in QDCSM3 is caused by the large nucleon size b used 
there. As b decreases, the AA attraction also decreases. 
This sensitivity to b is not seen in ChQMs, thus showing 
that it is caused by the QDCS mechanism. 

In fact, the 3 S^ N resonance disappears somewhere be- 
tween QDCSMI and QDCSM3. The critical value 6 crit 
below which the resonance disappears can be estimated 
under the assumption that the mass shift AM caused by 
the coupling to the NN channels is the same for all pa- 
rameter sets. The critical point then appears when the 
bare 3 S AA bound-state mass is 2186 MeV. Interpolation 
from the bound-state masses shown in Table IIIII gives 
6 crit w 0.53 fm. 

In contrast, the decrease of the pure 3 S AA bound state 



FIG. 2: (a) 3 Sf iv phase shifts calculated with two or five cou- 
pled color-singlet channels defined in Table ITTT1 in two different 
quark models ChQM2 and QDCSMI. (b) 3 Si N phase shifts 
calculated in the ChQM2 using different numbers of coupled 
channels, as described in more detail in the text. 



mass for ChQM2s as b increases from 0.52 fm to 0.60 fm 
is only 20% that of QDCSMs. Furthermore, b can be 
increased to only 0.645 fm, for the confinement potential 
strength a c turns negative above that value and hence no 
ChQM2s can be constructed. At b — 0.645 fm, the pure 
3 S AA bound state mass is 2332 MeV, which is too large 
for the system to resonate on coupling to NN channels. 
Hence ChQM2s have no 3 S^ N resonance. 

The pure 3 S AA bound state mass for ChQMl is 2274 
MeV, 92 MeV below the ChQM2 mass. If its b depen- 
dence is the same as ChQM2s, the decrease in bound 
state mass is insufficient, by about 50 MeV, to induce a 
resonance. 

Turning now to the NN phase shifts, Fig. shows 
that ChQM2 and QDCSMI give quite similar results, 
with the 5cc treatment giving much more attraction in 
both quark models, especially for E cm < 150 MeV. There 
is fair agreement with experiment for E cm < 150 MeV, 
but all theoretical phase shifts become increasingly too 
attractive at higher energies. 

Fig. H;b) shows the 14-channel 3 S^ N phase shifts 
calculated in ChQM2 with the addition of nine h.c. 
channels: eight 3 {S,D) 1 channels of 2 A 8 2 A 8 , 4 V 8 4 iV 8 , 
4 iV 8 2 V 8 , 2 V 8 2 7V 8 , and 7 L>i 4 A 8 4 iV 8 . Their inclusion 
causes the phase shift to become only a little more at- 
tractive. 

The figure also shows the phase shifts obtained af- 
ter the two arbitrary increases of the color confinement 
strength involving h.c. channels made previously for the 
3 D^ N system. The phase shifts are now noticeably dif- 
ferent from each other. Both are considerably more at- 
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TABLE IV: The dibaryon or resonance mass and decay width, 
all in MeV, in four quark models for 1 = 1 states. The chan- 
nels included in the J p = + state are lc ( 1 5 ( ^ A only), 2cc 
(lc + 1 S$' N ), and 4cc (2cc + 5 D^ A + 5 D$ A ). The channels 
included in the J p = 2 + state are lc ( 5 S2 A only), and 2cc 
(lc + 1 Z>f N ). The bound-state masses for ChQMl are 2304 
MeV for 1 Sq A and 2171 MeV for 5 5*f A 0. 





ChQM2 


ChQM2a 


QDCSMO QDCSM1 QDCSM3 




M 


M 


M 


M 


m r 




J p = 0+ 










lc 


2395 


2390 


2335 


2231 


2148 - 


2cc 


nr° 


nr 


nr 


nr 


2448 106 


4cc 


nr 


nr 


nr 


nr 


2433 128 




J p = 2+ 










lc 


ub 6 


ub 


ub 


ub 


2167 - 


2cc 


nr 


nr 


nr 


nr 


2168 4 



a No resonance in these coupled channels. 
'Unbound. 



tractive than those for case 14cc, but still not attractive 
enough at high energies for a resonance to appear below 
the AA threshold. These large changes in the low-energy 
phase shifts are inconsistent with experiment, thus show- 
ing that these additional h.c. effects can now be excluded. 

Two other NN partial waves merit a short discussion. 
The 5cc 3 D^ N phase shifts, like their 3 Si N partners, are 
nonresonant except for QDCSM3. All theoretical phase 
shifts agree well with experiment, to around 200 MeV. 
The quality of the theoretical 3 S"f N and 3 £»f N phase 
shifts shows that both quark models give good descrip- 
tions of the longer range part of the effective isoscalar 
central potential. The isoscalar 5 S AA is Pauli forbid- 
den, while the pure 3 D AA state is unbound. As a result, 
the 2cc 3 D2 N phase shifts are nonresonant. The QD- 
CSM values agree quite well with experiment, while the 
ChQMs agree less well. 



B. 1=1 states 

Table IIVI summarizes the results for two isovector 
states with possible resonances. The J p = + ( 1 S' ( ^ Ar ) 
state is qualitatively similar to the isoscalar J p = 1 + 
( 3 S^ N ) state since they are mostly different spin states 
of the same dibaryon pairs in the same relative orbital an- 
gular momenta. The pure 1 S^ A bound state in QDCSM3 
is pushed up from its unperturbed energy of 2148 MeV 
by 300 MeV on coupling to the 1 S' ( ^ rAr channel by the 
presence of a lower-mass state, the well-known slightly 
unbound x Sq N state. The perturbed mass is still small 
enough for the system to resonate below the AA thresh- 
old. 

In the remaining five quark models, the pure A A mass 
is 80 — 250 MeV higher. In each case, the strong coupling 
to the NN channel pushes the state well into the AA 
continuum, thus preventing a resonance from materializ- 
ing. Following the procedure used for 3 S AA , the critical 
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FIG. 3: NN 1 So phase shifts calculated with two or four 
coupled color-singlet channels in two different quark models 
ChQM2 and QDCSM1. 



nucleon size below which the resonance disappears in the 
QDCSM is found to be & crit w 0.56 fm. 

Figure [3] shows the non-resonant behavior of the 
coupled-channel x Sq N phase shifts for ChQM2 and QD- 
CSM1. These phase shifts become much more attractive 
than the experimental values from SP07 as the scatter- 
ing energy increases into the resonance region, the effect 
being more than twice as strong as the similar behavior 
in the 3 S^ N state. All the quark models studied here do 
not give enough short-range repulsion in these S-states. 

In the J p = 2 + state, the pure 5 S'^ rA bound state 
appears at roughly the same mass straddling the NA 
threshold in all six quark models used. (Model differ- 
ences in the NA 5-state masses are much smaller for 
the high intrinsic spin state than for the low intrinsic 
spin states, just like the model differences in the AA in- 
states.) The pure 5 S^ A mass is pushed up only a little 
by coupling to the 1 D2 N continuum. It remains bound 
only in QDCSM3. The pure 5 S 2 NA binding energy for 
ChQMl is only 0.14 MeV. So the state is unlikely to re- 
main below the iVA threshold after coupling to the NN 
channel. 

The pure NA state is unbound in ChQM2. (It be- 
comes bound if the attractive cental part of the scalar 
potential V shown in Eq. (J2J) is increased in strength by 
a multiplicative factor 1.7.) The calculated non-resonant 
NN phase shifts are shown in Fig. [4] for ChQM2 and 
QDCSM1. The prominent cusp is a threshold or Wigner 
cusp with its maximum located right at the NA thresh- 
old. The phase shift above the threshold is the phase of 
Six, where the subscript 1 denotes the NN channel. For 
comparison, the resonant phase shifts for QDCSM3 are 
also given. 

Even though a resonance appears only in one quark 
model in our limited theoretical treatment, the masses 
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FIG. 4: AW 1 Z?2 phase shifts calculated with two coupled 
color-singlet channels in three different quark models ChQM2, 
QDCSM1 and QDCSM3. 



involved are sufficiently close to one another and to the 
A A threshold so that they describe similar dynamical sit- 
uations to within the uncertainties of the models. More- 
over, the large A width when included would cause the 
state to straddle the A A threshold for all these quark 
models. We therefore consider a 5 S^ rA resonance near 
the A A threshold to be possible in all these quark mod- 
els. 

In fact, inelastic Argand looping (which we shall define 
in Sect. IIVI) has been obtained by Entem, Fernandez and 
Valcarce [23[ in the 1 D 2 and possibly also 3 F 3 systems 
for a ChQM having a s = 0.4977, only a little larger than 
the value 0.485 used in ChQMl or ChQM2. The crucial 
feature in their treatment is the explicit inclusion of A A 
inelasticities by giving decay widths to the As appear- 
ing in the coupled-channel treatment. There exist quite 
extensive PW solutions of both A A [l6| and nd [24[ scat- 
tering amplitudes in these and other isovector dibaryon 
systems. They could yield interesting information con- 
cerning quark dynamics in this resonance region. 

We do not find any resonance attributable to an A A or 
A A bound state in any of the quark models in the follow- 
ing four isovector states: (a) the 3 Pq i, 3 F3 states (with 
each state calculated using three coupled color-singlet 
channels of the same quantum numbers for NN, NA 
and AA constituents, respectively), and (b) the J p = 2~ 
state (using the four color-singlet channels 3 P 2 of NN, 
NA, AA, and 7 P 2 AA ). 

It is worthwhile to show in Fig. [5] the difference be- 
tween ChQM2 and QDCSM1 in the well-known decom- 
position into central, spin-orbit and tensor components 
phase shifts: 



of the 3 Pf N 
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We see that the inclusion of one pion exchange (OPE) 
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FIG. 5: The central, spin-orbit and tensor components of 
the 3 pj VJV phase shifts calculated with three coupled color- 
singlet channels (containing NN, NA and AA) in two differ- 
ent quark models ChQM2 and QDCSM1. 



takes good care of the tensor component, but both quark 
models give too weak spin-orbit components and too 
repulsive central components, especially for QDCSM1 

ai.. 

As is well known [2j, the problem with the spin-orbit 
component comes about because the included OPE po- 
tential though clearly needed to generate a AA tensor 
force also contributes to the A— A mass difference. The 
color coupling constant a s needed to account for this 
mass difference is then reduced to 0.3-1.0 from the value 
1.7 in quark models without pion exchange. This weaker 
a s gives in turn a weaker baryon-baryon spin-orbit po- 
tential from OGE. The additional spin-orbit contribu- 
tion from the scalar exchange used in the ChQMs is not 
enough to compensate for the deficit. Resolution of this 
problem would require significant modifications to the 
quark models used. 

A simple way to study the spin-orbit problem in the 
present limited objective of looking for dibaryon reso- 
nances without overhauling these quark models is to just 
modify a term or add terms to the quark-quark inter- 
action phcnomcnologically. A number of related quark 
models are thus generated: a modified ChQM2m model, 
where the one gluon exchange strength a s of the spin- 



orbit potential V, 



G,LS 



has been increased 5 times, and 
a number of ChQM2p(f) models with the additional p 
meson exchange potential V displayed in Eq. ([T]). Here 



/ = 



&chv0 



(8) 



is the multiplicative increase of the p-quark-quark cou- 
pling strength above the usual and customary value 
o-chvo — 0.021, corresponding to the coupling constant 
9chv0 = 2.351 used in [T3|. Since the effects on these 
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FIG. 6: The central, spin-orbit and tensor components of the 
3pNN pjjg^gg shifts calculated with only a single uncoupled 
NN channel in the quark model ChQM2p(ll). 



P-wave NN phase shifts of the coupling to NA and 
AA channels are quite small, we use only a single color- 
singlet NN channel to calculate the NN phase shifts for 
these modified models. Fig. [5] shows that the resulting 
central, spin-orbit and tensor components of the 3 p^ N 
phase shifts for the quark model ChQM2p(ll) give quite 
good agreement with the experimental SP07 values. By 
diagonalizing the Hamiltonian matrix in the appropriate 
NA channel, we find that none of the modified quark 
models described in this paragraph has a pure NA P- 
wave bound state. 

The difficulty of forming P-wave bound states can 
readily be appreciated in the attractive square-well po- 
tential model. It is well known that to bind a P-state, its 
attractive potential dep th must be four times that needed 
to bind an S-state [25|]. To see how far we are from NA 
P-wave bound states for our quark models, we increase 
the strength of the attractive central scalar potential V 7 
in ChQM2 by a multiplicative factor f s . The 3 P 2 JVA state 
then becomes bound at f s > 3.1, with a binding energy 
of 6 (22) MeV when f s = 3.2 (3.4). The potentials in the 
other 3 pj^ A states are less attractive. 

A similar situation holds for ChQM2p(ll): Its 3 P 2 JVA 
( 5 P 2 WA ) state is unbound. It becomes bound only when 
its attractive central scalar potential strength has been 
increased by a multiplicative factor f s w 3.0 (2.6). It 
is thus clear that NA P-wave bound states would ap- 
pear only for quark models with substantially stronger 
attraction than those studied here. 

Concerning the missing P-wave attraction in both 
ChQM2 and QDCSM1, one cannot just use the isoscalar 
scalar S meson of mass 980 MeV that appears in the Bonn 
potentials [26| because of the cutoff mass A « 830 MeV 
used in our chiral quark models. It is not a trivial prob- 
lem to reconcile these two classes of models for nuclear 



Turning now to our experimental knowledge of pos- 
sible resonances in NN P-statcs, wc recall that in the 
PW analysis FA91 [HI], resonance poles are found for the 
isovector odd-parity NN states 3 p2, 3 P2 and 3 p3. These 
empirical resonance-like solutions reproduce the empiri- 
cal Argand loopings of the PW solutions, but many stud- 
ies in the past [27[ have left unresolved the question of 
whether these Argand loopings represent new dibaryon 
resonances. The difficulty centers around the observation 
of Brayshaw [28[ that when decay channels with three or 
more final particles are present, Argand looping can ap- 
pear in models known to have no resonance because their 
S-matrix has no pole. Brayshaw has given an explicit dy- 
namical example for the case of the pp 1 Di state coupled 
to NA and ir + d channels. The strong energy dependence 
that causes the Argand looping in his model comes from 
a logarithmic rescattering singularity in the iVA «-> n + d 
transition amplitude near the NA threshold. The phys- 
ical situation this singularity describes is the oscillation 
or exchange of a nucleon between the decaying A and 
the second or spectator nucleon with which it forms the 
bound deuteron d. In Brayshaw's model, there is no new 
1 D 1 2 P resonance near the iVA threshold. 

To complete our very brief review of resonance condi- 
tions, we should mention that recent studies of irN reso- 
nances [29j, |30j have shown that the speed test can track 
the positions of resonance poles, if present, more reliably 
than the time delay criterion. (The speed test determines 
the resonance energy from the maximum speed of Argand 
looping as a function of the on-shcll kinetic energy, while 
the time delay test locates it by maximizing the positive 
time delay of the scattered wave packet relative to the 
free wave packet.) A different kind of complication can 
appear in dynamical models that already have resonance 
poles, namely that the speed test can fail because the 
Argand looping does not have a solution with maximum 
speed 29]. This is an extension to another dynamical 
model (one having two coupled nonrelativistic two-body 
channels) of another old observation of Brayshaw that 
when relativistic many-body channels are present, true 
resonance poles can appear without any Argand looping 

S2. 



In view of all these complications, we shall take the 
tentative but conservative position that a promising 
dibaryon resonance is one involving at least one A or N* 
baryon that is a bound state below the dibaryon breakup 
threshold in the absence of a centrifugal potential when 
these baryons are treated as stable particles. In prac- 
tice, the only excited baryon we are able to describe with 
some degree of confidence is the A. In the limited context 
of our quark models, we consider the 1 D^ !N structure a 
promising dibaryon resonance, at least for some of our 
quark models, but not the NN P-wave Argand loopings. 
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TABLE V: Mass, decay widths (both in MeV) and branching 
ratio of theoretical baryon resonances in five quark models and 
comparison with partial-wave analyses of experimental data 
from SP07 [l]J and FA9l[H. Mr for ChQM2a is estimated 
from the 4cc value and is shown within parentheses. 
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3 Si 
3 D 1 


L D 2 


■ i D 3 
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Si 
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AT A 
iV A 
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Do 










ChQM2: 








Mr 








2393 


YnN 








14 


Tinel 








136 


Bnn 








0.09 


Mr 


ChQM2a: 






(2404) 
9 


Tinel 








149 


Bnn 








0.06 




QDCSMO: 








Mr 








2400 


r jviv 








14 


Tinci 








144 


Bnn 








0.09 




QDCSM1: 








Mr 








2357 










14 


Tinel 








96 


Bnn 








0.13 




0.60 


0.57 








QDCSM3: 








Mr 


2433 


2393 


2168 


2273 


Ynn 
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70 


4 


17 


Tinel 


190 


136 


117 


33 


Bnn 


0.41 


0.34 


0.03 


0.34 




SP07: 








Mr 


none 


none? 


2148 a 


? 


r 






118 a 




Bnn 






0.29 


0.26 6 



"Pole position of FA91. 
6 At W = 2400 MeV. 



IV. DISCUSSION 

The theoretical resonance properties calculated in the 
last section are summarized in Table [V] The first line 
for each dibaryon type gives the dominant PW. The PW 
responsible for the resonance trapping in the theory is 
shown in bold type. The experimental information used 
for the comparison is the PW solution SP07 [lj| of NN 
scattering data. The four states are arranged in order of 
increasing relative orbital angular momentum £nn- 

Each of our calculated resonances is an elastic res- 
onance where the scattering phase shift rises sharply 
through 7r/2. The Argand plot of its complex PW am- 
plitude 



T 



5-1 

2i 



(9) 



shows rapid counterclockwise motion on the unitarity 
circle. This mathematical behavior describes a physi- 
cal picture where the NN system resonates or continues 
to "sound" due to its partial trapping into the resonance 
region, namely the closed channel containing one or two 
As. Each of these elastic resonances has a finite but very 
small elastic (or NN) width. 

The resonance properties shown in the table must be 
corrected for the width of a decaying A, leading to the 
appearance of pionic channels. Inelasticities cause the 
reduction \S\ < 1 and the restriction of the Argand plot 
to the interior of the unitarity circle. We shall avoid 
Brayshaw's two complications from many-body channels 
by considering only channels where for stable As, the 
dibaryon system in our quark model treatment is a bound 
or almost bound two-body state. This restriction al- 
lows us to take the standard position that for these spe- 
cial states, rapid counterclockwise Argand looping is an 
acceptable signal of an inelastic resonance [27|, |32|. In 
physical terms, leakage of the trapped system into pionic 
channels reduces the effect "heard" in the NN channel 
but does not eliminate it altogether. From the perspec- 
tive of the quark models used here, the I — 1 odd-parity 
NN Argand structures are not promising candidates for 
dibaryon resonances. 

We shall estimate inelasticities only at the crudest level 
of branching ratios, with 



B 



NN 



NN 



(10) 



where Y 



NN 



depends on the inelastic width 



Tinci caused by decaying As. Close to the breakup 
threshold where the As are almost on-shell, the in- 
elastic width can be related approximately to the A 
width Yt A = 120 MeV in free space by only accounting 
for the reduction in phase space available to a decay- 
ing bound A whose mass has been reduced to roughly 
M&a ~ MfA — Q.5.B. where B is the binding energy of 
the dibaryon. Then [33j 



Y bA (M bA ) « r 



/A 



k 2 b e p(M bA ) 
kfp(M fA y 



(11) 



where k is the pion momentum in the rest frame of the 
decaying A, £ = 1 is the pion angular momentum, and 



p(M) 



kE^EN 
M 



(12) 



is the two-body decay phase space at mass M when each 
decay product has cm. energy E i: i = ir, N. We shall use 
this crude estimate indiscriminately even far below the 
breakup threshold, but the harm done is not great be- 
cause most promising resonances are near the threshold. 

If each decaying A in the dibaryon has a Breit-Wigner 
(BW) distribution of width Yb A , the total mass of two 
decaying As can be shown to have a BW distribution 
with width 2Yb A . Hence 



Tmci ~ n b ArbA(A4 A ), 



(13) 
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ribA is the number of As in the resonance. The results for 
different theoretical dibaryons are summarized in Table 
IVl The shift in resonance masses caused by the coupling 
to pionic channels has not been included in these esti- 
mates. 

The experimental branching ratios shown in Table [V] 
are obtained from the energy-dependent SP07 PW solu- 
tion [lfij at the stated resonance masses using the formula 
for coupled PWs 

d _ <7«u _ |Tu| 2 + |T 12 | 2 

DNN - - f—^ , (14) 

ftot,J 1ml n 

while T\2 = for an uncoupled PW. Here the subscripts 
i,j in Tij are channel labels. 

Our theoretical estimates can now be compared to the 
experimental results from the PW analysis of NN scat- 
tering data, in Table |V] The absence in the SP07 PW 
solution of a resonance accessible from the 1 Sq FN channel 
and its probable absence in the 3 S^ N give an approxi- 
mate upper bound on the nucleon size of 6 crlt « 0.53 fm 
for the QDCSM. This upper bound causes the 3 D^ N res- 
onance mass to exceed 2340 McV, the interpolated value 
for the QDCSM at b = 0.53 fm. 

The 3 D^ N resonance appears also in both ChQM2 
and ChQM2a at around 2400 MeV, the resonance mass 
being not very sensitive to b. The resonance mass has 
not been calculated for the ChQMl, but is about 2420 
MeV. NN S'-wave resonances do not appear in ChQM 
with quadratic or linear confinement. 

Table fVl also shows that the theoretical 3 D^ N branch- 
ing ratio Bnn is smaller than the SP07 value. This 
means that the theoretical coupling to the trapping AA 
channel is too weak. The theoretical branching ratio 
Bnn for the QDCSM3 resonance in the 1 D^ N state is 
much smaller than the SP07 value except for QDCSM3. 
This suggests that the theoretical coupling to the trap- 
ping NA channel is also too weak. 

A. The ABC effect 

The ABC effect, named after Abashian, Booth and 
Crowe 34] who first observed it, describes an enhance- 
ment above phase space of the missing-mass spectrum 
of the inelastic fusion reaction pd^ 3 HcX . Subsequent 
experimental studies have been reviewed recently by 
Clement et al. (or WASA07) [11], who also report prelim- 
inary results for the exclusive reaction pi— > 3 He7T7r from 
the WAS A Collaboration, and by Bashkanov [35| . 

The enhancement is associated with a tttt invariant 
mass < 340 MeV, with the two pions emitted in parallel 
and in relative S'-wave opposite in direction to the recoil- 
ing nucleus. The structure is isoscalar because it is seen 
in dd^ 4 HeA°, but not in dp^ 3 YLX+ . 

In the reaction pd— > 3 He7r°7r°, the differential cross sec- 
tion da/dM w o„o at fixed M3, the invariant 3 He7r° mass, 
has a maximum at w 3080 MeV, just under the value of 
2M N + M A « 3100 MeV. The cross section has a FWHM 



TABLE VI: Experimental peak position (in MeV) and pro- 
duction cross section admr (in mb) and width (in MeV) of the 
ABC effect compared to the estimate from the partial-wave 
solution SP07 of np scattering amplitudes under the assump- 
tion that the ABC effect comes from a dibaryon resonance 
in the specified PW at the resonance mass 2410 MeV. These 
estimates are described in more detail in the text. 



Reaction 

Ref. 

PW 


np — > dX np — > dniv 
H-TA73[36] WASA07[11J 


np — > np 
SP07[16] 
3 S"i 3 D 3 


Peak 


2460° 


2410 


none? ? 


Ctot, J 






3.6 3.2 


Cinel.J 






0.41 2.3 


&dnn . J 


0.63 


2.3 


wO.12 R!0.8 


r 


> 220" 


< 100 


160 6 160* 



"Theoretical calculation of [37| . 
'From Eq. (TTJ at 2410 MeV. 



of about r 3 ps 130 MeV, the same as T/a- (These num- 
bers are from Fig. 5-5 of [Hj].) Hence the pion in M3 
appears to have come from the decay of a slightly bound 
A. These features of the experimental data suggest that 
the ABC effect in these reactions comes from the decay 
of a AA bound state. The preliminary WASA07 res- 
onance peaks at 2410 MeV, about 50 MeV below the 
AA threshold, with a width < 100 MeV. If the estimated 
width holds up, it would eliminate the larger values of 
2r bA « 160 MeVffrom Table EQ) to 2r 3 w 260 MeV 
(from Fig. 5-5 of (3a|). As pointed out by WASA07, such 
an outcome would disagree with the situation in the D2 
resonance whose width is close to the free A width even 
though the resonance straddles the NA breakup thresh- 
old. 

The energy dependence of the production cross sec- 
tion Ud-rnr (all final pion states) has been measured. Two 
rough fits to different data are shown in Table IVT1 The re- 
cent preliminary WASA07 results [ll| agree rou ghly with 
the older Heidelberg-Tel- Aviv data (H-TA73) [36(. In- 
formation can also be deduced from NN scattering. For 
comparison, the table gives the PW total and inelastic (or 
reaction) cross sections from the latest energy-dependent 
solution SP07 [lq ]. We must next estimate the fraction 
of (Tinei that goes through the dirir channel. 

The estimate is made by using the following two as- 
sumptions: (1) The single-pion production cross sec- 
tion in either J state is certainly not zero because a 
pion can be produced with the dibaryon left in isovec- 
tor states. However, in order to maximize our estimate 
for Ud-K-n, we shall ignore the contributions of all one- 
pion decay branches. (2) The experimental cross section 
<7 ( j 7r + 7r - (<J n pn+TT-) is known to be 0.270 (0.55) mb at 
W = 2340 MeV, and 0.33 (4.05) mb at W = 2510 MeV 
H]. Thesepoints straddle the 2410 MeV ABC peak of 
WASA07 [11(. Both cross sections should increase sig- 
nificantly as one approaches the ABC peak from below. 
Above the ABC peak, the larger dibaryon breakup cross 
section probably reflects an increase in phase space, in- 
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eluding the number of other contributing states. We now 
assume that both cross sections at and below the the 
ABC peak are dominated by the same resonance in one of 
the two J states in the table. Using the ratio 0.27/0.55 ~ 
0.5 of these cross sections at W — 2340 MeV for the en- 
tire resonance, we therefore assume that a dun ~ 0mei/3 
to get the rough and perhaps generous estimates shown 
in the table. Finally the estimated width for the two AA 
states is that of Eq. ([13]) which includes the reduction in 
phase space for bound As. 

Table [VTl suggests that it is relatively unlikely that the 
ABC effect originates from a dibaryon resonance in the 
3 S^ N channel. The main reason is that the NN scat- 
tering described by SP07 is highly elastic so that u,i„ 
is far too small. The situation for the 3 D^ N channel 
is more promising but not without difficulty unless the 
preliminary WASA07 estimate of is reduced. 

Additional information can be obtained from the en- 
ergy dependence of the SP07 3 D^ N scattering amplitude 
[16l | . They are non-resonant at the QDCSM3 resonance 
mass of 2270 MeV, but they are too uncertain at the 
ABC peak at 2410 MeV to settle the question of an NN 
resonance there. 

The resonance widths are also of interest, especially for 
the exclusive reaction that excludes contributions from 
three and more pions on the high-energy side of the pos- 
sible resonance. The theoretical calculation of Bar-Nir, 
Risser and Schuster [37[ is based on a one-pion exchange 
excitation to two As followed by a pion emission from 
each A. Their calculated resonance width is close to free- 
space value 2TfA = 228 MeV used in their calculation. 
Our estimated decay width shown in Table I VII is much 
smaller but not as small as the preliminary WASA07 
value. As for the branching ratio Bnn given in Ta- 
ble El the calculated value for our quark models seems 
too small, but the experimental value from SP07 at 2400 
MeV is not necessarily reliable. 

In his study of the ABC effect, Alvarez-Ruso [39| 
has pointed out that the AA contribution is greatly re- 
duced when short range repulsive correlations are in- 
cluded in the NN channels. Then the cross section at 
W = 2240 MeV, some distance below the ABC peak, is 
found to be dominated by the AW*(1440), with both 
pions emitted by the decaying Roper resonance. How- 
ever, at this lower energy region, the SP07 np 3 D$ Ar- 
gand phase motion is not resonant. 

Short range correlations are already included in the 
quark models used here. They are not the short-range re- 
pulsion from the exchange of vector mesons (specifically 
the isoscalar u> meson) , which would reduce if not elim- 
inate the AA resonance. Our short range correlations 
come from Pauli antisymmetrization and channel cou- 
pling effects generated by overlapping clusters of quarks, 
including baryon excitations and hidden-color configu- 
rations that enhance rather than reduce these short- 
distance phenomena. 

The relative importance of these explicit quark effects 
found in our studies also comes from the use of large 



baryon clusters of quarks in all the quark models used 
here. The situation could be different if the baryon 
"bags" are small and the meson clouds around them 
are thicker. Our calculated results for the dibaryon spec- 
trum is sensitive to the model nucleon size used in the 
QDCSM, but not in the ChQM. The baryon spectrum 
on the other hand is quite sensitive to the model nucleon 
size, especially in the radial excitations. Past calcula- 
tions in the ChQM Q favors the choice near b = 0.52 fm. 
With this choice, the theoretical 3 D^ N resonance arising 
from the 7 S^ A bound state appears at about 2390 MeV 
in ChQM with quadratic confinement, probably at 2420 
MeV with linear confinement, and at 2360 MeV for the 
QDCSM. 

Experimental ed form factors at large momentum 
transfers that show the premature dominance of six- 
quark effects seem to favor the kind of quark models 
studied here over the more traditional short-range NN 
repulsive correlation [4l| . The experimental confirmation 
of a NN 3 D% resonance would be a dramatic demonstra- 
tion of quark effects in the resonance region. Its experi- 
mental non confirmation on the other hand would point 
to a missing short-range repulsion in our quark models, 
a repulsion that is usually attributed to vector meson ex- 
changes in traditional meson exchange models of nuclear 
forces. 



V. CONCLUSION 

We have studied resonances in NN scattering in a the- 
oretical treatment of two baryon clusters of quarks in- 
teracting by Pauli antisymmetrization and by gluon and 
pion exchanges. The nucleons can resonate by changing 
into As, but only if the resulting baryons attract each 
other with sufficient strength to stay below its 5-wave 
breakup threshold. The absence of NN S-w&ve reso- 
nances in the SP07 partial-wave amplitudes places an 
approximate upper bound of b < 0.53 fm on the nucleon 
size in the QDCS quark model. This restriction in turn 
requires that the 3 D^ N + 7 S AA resonance mass should 
exceed 2340 MeV. In ChQMs, the NN system does not 
resonate in relative-5 waves, but it has a 3 D% resonance 
at 2390-2420 MeV. This 3 D 3 resonance is thus a promis- 
ing candidate for the explanation of the ABC structure at 
2410 MeV in the production cross section of the reaction 
pn— >G?7T7r. 

The most promising isovector NN resonance candi- 
date in our quark models appears in the 1 D^ N state and 
comes from a bound or almost bound 5 S AA state. None 
of the quark models used has bound VA P-statcs that 
might generate odd-parity isovector resonances. 

It is satisfying that these simple quark models con- 
taining only a few adjustable parameters fitting the N, A 
masses and the deuteron binding energy can yield phys- 
ically interesting information about the possibility of 
dibaryon resonances at the much higher energies near the 
AA and NA thresholds of NN scattering. Their success 
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is partly due to the fact a good part of the available colli- 
sion energy in the center of mass frame has been used to 
excite the nucleons into one or two As whose mass is fit- 
ted by the models. The system is thus effectively at much 
lower energies in these A channels. The success is also 
derived from the ability of these simple models to capture 
some essential features of the baryon-baryon interactions 
in different energy regimes in the many channels involved 
in the calculation. 

The quark models used have many shortcomings. In 
the context of nonrclativistic models alone, a quantita- 
tive fit to NN phase shifts appears difficult without the 
fine tuning provided by the addition of the many me- 
son exchange terms that appear in conventional boson- 
exchange potentials The need to use both quarks 
and meson exchanges suggests that the resonance region 
of interest in this paper is a transition region between the 
traditional low-energy regime of baryons interacting via 
meson exchanges and the high-energy regime of quarks 
interacting by quantum chromodynamics. 



From a more technical perspective, our theoretical de- 
scription can be improved by treating the As as decaying 
particles. It would be difficult to go beyond this im- 
provement because explicit pion channels contain three 
or four bodies. In our discussion of the ABC effect, it is 
of considerable interest to improve upon our very rough 
estimate of the partial-wave pn^dirir production cross 
sections from NN partial-wave amplitudes. In spite of 
these limitations, it is clear that additional experimental 
knowledge and theoretical studies of NN properties in 
the NN resonance region near the NA and AA thresh- 
olds will add significantly to our understanding of quark 
dynamics between baryons. 

Finally we should add that the final report of the 
CELSIUS- WASA Collaboration on the ABC effect has 
now appeared [43| . The structure in the total cross sec- 
tion for the pn — > ci7r + 7r~ reaction centers at 2.39 GeV 
with a width of 90 MeV. 

This work is supported by the NSFC grant 10375030, 
90503011, 10435080. 
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